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Precis 


The application of boundary layer suction is as a rule con- 
cerned with single slits, while the majority of the theories 
deal with the case of continuous suction. In this report, a 
simple theory for suction by single slits is developed. The 
theory is connected with the calculation procedure for lami- 
nar and turbulent boundary layers, described in DVL Reports 
8i;(1959) and 136(1960). 

In this theory there are only two unknowns for the velo- 
city profile, a thickness parameter and a shape parameter. The 
basic concept of suction theory now consists in ascertaining the 
values of these tv;o parameters behind the suction slit from the 
parameter values in front of the slit, since the arriving velo- 
city profile with the given suction quantity is "clipped" ac- 
cording to the continuity condition and the parameter values for 
the remaining profile are specified in an appropriate manner. 

The calculation then continues behind the slit with this para- 
meter value. 

The theory is fully developed in this report at first only 
for turbulent boundary layers on the basis of velocity profiles 
according to the common and simple exponential law. As early as 
1943 the author had presented an incomplete form of this theory 
in a limited-access report of the Aerodynamical Research Insti- 
tute, Go'ttingen, and compared it with measurements. Several 
essential findings of this older report have been reproduced in 
the current report. 
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APPROXIMATION THEORY FOR BOUNDARY LAYER SUCTION 
THROUGH INDIVIDUAL SLITS 

A. ¥alz 


I, Introduction . 

The viscosity of fluids has effects that are almost always 
undesirable in flow technology. Among these are, e.g., the 
resistance of bodies in the flow to friction and pressure, the 
pressure losses in pipelines, and the detachment of the flow 
with its subsequent disadvantages. For the order of Reynold's 
numbers prevalent in flow technology, the viscosity is practi- 
cally restricted in its action to a very thin boundary layer. 

Thus the above unfavorable phenomena are also to be referred to 
processes in the boundary layer, for large Reynold's numbers. 

L. PRANDTL Z”i7 first to show (190ij.) the connection 

between the processes in the boundary layer and the outer flow 
of potential. From then onwards, the way was clear for all in- 
vestigations that attempted to choose the contours of objects in 
a flow or the cross sectional curve’ of channels in such a way 
that the boundary layer processes, in respect of flow resis- 
tance and detachment phenomena, occur in an optimal manner. Some 
of this work led to the development of airfoil profiles with 
small resistance (the laminar profile) and at the same time a 
high maximxim lift Recently, R, EPPLER /,”*67 has 

treated this problem with the modern methods of potential the- 
02 ?y and boundary layer theory. F. X. Wortman has made ex- 
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tensive measurements in regard to this problem. But also in 
this connection we should mention all the mathematically profound 
and experimentally difficult works that explore the conditions 
for stability of the laminar boundary layer (cf. the newest 
findings of H. GSrtler’s research group, Freiburg Institute, 

V/est Germany, and for a bibliography see, e.g,, DVL Report 8), 

Let us now assume that we have exhausted all possibili- 
ties to constitute an optimal pressure, working on the boun- 
dary layer from without. It is now to be asked whether there 
are yet other physical or technical possibilities for inter- 
vention in the boundary layer processes. In critically evalu- 
ating these possibilities, it is seen that a distinction should 
be made between the measures for preventing flow detachment 
and those for reducing the frictional resistance . This will 
be illustrated with the help of figures 1 through ?, 


II, The Fundamental Possibilities for Influencing the Boundary ^ 
Layer , 


II, 1 , Measures to Prevent Stream Detachment , 

The manner in which a boundary layer (assumed to be lami- 
nar) affects a pressure gradient dp/dx in the flow direction 
should be immediately apparent from fig, 1, In the case of 
dp/dx = 0 (flow along a level plate), we find something like 
the indicated curve of the velocity u with the distance y from 
the wall (Blasius profile /f^Z). The particles of fluid close 
to the wall still have a slight kinetic energy, by which they 
become "playthings” of external pressure gradients dp/dx, which 
are nonetheless active or "impressed” in the neighborhood of the 
wall. If dp/dx ■< 0, then the particles near the vjall are ac- 
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celerated. Velocity profiles that are more full than the Bla- 
sius profile are the result. In the case of a positive pres- 
sure gradient (dp/dx>0), kinetic energy is expended in the 
formation of the static pressure field. The velocity close to 
the wall decreases even further than for the Blasius profile. 

The result is velocity profiles with a turning point so that, 
in the neighborhood of the wall and for sufficiently large dp/ 
dx, there may occur an inversion of the velocity direction, with 
consequent flow detachment. 

To avoid flow detachment in this case we can use, for ex- 
ample, a wall stage with adjacent suction channel, as shown 
schematically in fig, 2, The "sick" flow of the boundary layer 
near the wall disappears now in the suction channel (naturally 
one must provide for a lower pressure in the channel with res- 
pect to the outer pressure), and a new, "healthy" boundary lay- 
er flov; begins behind the suction point. It has a fuller ve- 
locity profile, able to overcome further pressure rises. 

A comparison of the velocity profiles before and after the 
suction point shows immediately that, while the suction has pre- 
vented detachment of the flow, it leads to a greater wall fric- 
tion behind the suction point, due to the greater velocity. V7e 
thus cannot expect to reduce the frictional resistance by suc- 
tion alone. The points later discussed in connection with figs, 

8 and 9 have some bearing here, V/e now consider still other 
possibilities for preventing flow detachment. 

It is technologically more simple to produce a suction by 
a slit in the surface, as sketched in fig, 3, The result is 
practically the same as the case in fig. 2; the energy-poor /7 
boundary-layer flow beneath the dot-and-dash streamline is en- 
gulfed in the suction slit. Behind this slit, a fuller and en- 
ergy-rich velocity profile is again available to confront new 
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pressure increases vrithout danger of detachment. 

Fig. i|. represents the case of a continuous suction, such 
as by a porous wall. It is the object of much research (see 
the bibliography in shall content ourselves here 

with a reference to the bibliography. Another possible way 
to increase the kinetic energy of the boundary layer near the 
wall is by blowing , tangential to the surface and in the flow 
direction, as shown in fig. By this means, the velocity 
profile behind the blowing point becomes more full. But at 
the same time, the thickness of the boundary layer can grow, 
which is not the case for the preceding measures. 

Finally, it is possible to influence the boundary layer 
in prevention of flow detachment by moving the wall in the flow 
direction. But this measure has till now hardly been applied 
in practice, apparently due to the mechanical difficulties, 
(Example: the Flettner Rotor for producing a circulation flow, 

see /~107, p. 2I4.7). 

To illustrate what advantages can result to flow technol- 
ogy By the prevention of flow detachment, the classical experi- 
ment of L. PRANDTL on a short diffuser with two suction 

slits should be recalled, figs. 6 and 7. Also the successful 
attempts of the AVA Gottingen to increase the maximum lift of 
airfoil profiles should be consulted (for a summary of which 
see, e.g., 2fl27. 


11,2. Measures to Reduce the Frictional- Resistance (Maintain - 
tenance of the Boundary Layer* s Laminarity ) , 

For the measures depicted in figs, 1-5 we can say in sum- 
mary that, near the wall, the velocity within the boundary la- 
yer is increased. In this way, the flow is prevented from de- 
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tachraent in regions of pressure rise. But these measures, as 
already indicated, at first appear to be ill-suited to reduc- /8 
ing the frictional resistance, since the shearing stress of 
the wall increases, after Newton’s theorem, with the velocity 
gradients perpendicular to the wall (y-direction) . 




w 


dy /y=o 


( 1 ) 


u = the velocity component parallel to the wall (x-direction) 
y = the viscosity coefficient 


But reducing the thickness of the boundary layer and making the 
velocity profile more full can have an additional important ef- 
fect, about vjhich we have not yet spoken: the laminar condition 

of the boundary layer can be stabilized or — in other words — can 
be prevented from entering :he turbulent condition at turning 
points. The consequence of this is illustrated in figs. 8a, 8b, 
and 9. If the laminar and :he turbulent velocity profile are 
compared, figs. 8a and 8b, with the thickness of the boundary 
layer being equal in both cases, 6] “ > and the velocity at 

the edge of the boundary layer also being equal for the two, it 
is seen that the velocity gradient (3u/3y)y»o of a turbulent 
boundary layer is much grea-ter than that of a laminar, since the 
turbulent exchange movements bring kinetic energy into the wall 
proximity.^ In this way thk frictional resistance can be much 
greater for a turbulent, than a laminar boundary layer. If we 
consider, for example, the total frictional coefficient Cp of 
a level plate (W = total resistance per unit of width, 1 = length 
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If it is intended to prevent flow detachment, then the conver- 


sion of the boundary layer 
condition is to be regarded 
have as favorable an effect 


rrom the laminar to the turbulent 
as a desirable process, which can 
as suction. Thus, for flow proces- 


ses with a small Reynold's number, artificial turbulence, such 


as produced by a trip wire. 


tachment, often used in practice. 


is a simple measure to prevent de- 
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of the plate, J^= density of the fluid) 



W 



2 


) 



tw 


/ 

dlj-)=2jcfd(f) 

O 


with 





= the local coefficient of friction. 


then we find Cp and the Reynold* s number 





connected by the curve Cp(R^), which is very different for a 
laminar and a turbulent boundary layer (fig. 9). When the Rey- 
nold's numbers are greater than 10', (Cp)-turb almost one 

power of ten greater than 10*^, the differ- 

ence between the two can become greater still. 


If special measures are not taken to maintain the laminar- 
ity of the boundary layer, there occurs a transition at about 

(x = path length) from the laminar to the turbulent 

A 

condition, for flow along a smooth and even plate. This is a 
consequence of instability of the laminar boundary layer. Now 
it is stated in the stability theory of a laminar boundai?y la- 
yer that the layer can remain laminar, even at R > 10^, if 

a) the thickness of the boundary layer is reduced by suc- 
tion, 

b) the velocity profile is made more full than the Bla- 
sius profile of the level plate. 

However both of these are conditions that can be fulfilled by 
suction of the boundary layer. Apparently it is very impor- 
tant not to suck off more fluid than is absolutely necessary 
to maintain the stability of the laminar layer, since both mea- 



sures, a) and b), increase the wall shearing stress of the la- 
jjjj_jiar boundary layer* Furthermore, a certain amount of power 
must be expended in transporting the suctioned material, which 
increases with the quantity transported. Since the maintenance 
of boundary layer laminar! ty can improve the frictional coef- 
ficient by several orders of magnitude (especially for the case 
of large Reynold's numbers Rj^), these unfavorable consequences 
of suction play only a minor role in comparison to the theo- 
retically minimal quantity of suction. It is also apparent that 
the theoretically minimal quantity of suction is attained for 
the case of continuous suction (fig. i/.), which is also a mea- 
sure for preventing flow detachment. The technical realiza- 
tion of this case is naturally very difficult, since porous 
surfaces are easily occluded. 

As H, HOLSTEIN was the first to show, and as W, 

PPENNINGER has recently demonstrated with measurements, 

the optimal case of continuous suction can be approximated with 
a rather large number of slits, succeeding each other in the 
flow direction. The technique of suction through single slits 
is also promising for aircraft design, as seen from both past /1 0 
and recent experiences. 

, W. PFEHNINGER has briefly reported on new experimen- 

tal results of suction by single slits, during which it was pos- 
sible to maintain the laminar boundary layer on airfoils up to 
very high Reynold's numbers, 3 The coefficient of 

profile resistance €„ (with the power expenditure for the sue- 
tion taken into account) was Cp<^ 10 (For a turbulent boun- 
dary layer, Cp was nearly 6 times as great). The results, dis- 
cussed in ^~167, as to the amount of suction that should be em- 
ployed, are also in good accordance with the theoretical cal- 


2 

Gf, the bibliography on the tests in the AVA Gfittingen and the 
ETH Zurich, /“127 to 
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culations for the stability problem. 

According to a report by H, R, Head results have 

lately been achieved in England that are equally favorable. 

Thus we can consider the maintenance of boundary layer la- 
minarity by suction to be a special modern problem in flow me- 
chanics. 


III. Approximation Theory for Suction through a Slit . 

III.1 . A General Statement of the Problem . 

It is a great deal more easy to implement suction with sin- 
gle slits, than continuous suction through porous materials (dan- 
ger of occlusion). Yet for the theoretical treatment of this 
problem, we find a totally different picture, since we can easily 
specify exact solutions of the boundary layer equations, when 
dealing with the continuous suction of a laminar boundary layer 
(cf., e.g., B. H. Schlichting ^”"l07); but the theoretical treat- 
ment of single-slit suction is complicated (cf., e.g,, W. Wuest 
/“18, 197, W. Rheinboldt /“207). 

In this report we shall describe a simple theory for suc- 
tion with single slits. It can be applied to laminar, as well 
as turbulent boundary layers. This theory follows very closely 
the approximation theory that was presented in DVL Reports 61 ^ 
and 136 . The theory enables the design engineer to calculate, 
with sufficient accuracy, the amount and power of suction needed, 
for example, to raise the maximum lift coefficient of a given 
airfoil profile by a certain amount, or to prevent flow detach- 
ment in a given short diffuser. In applying this theory to la- 
minar boundary layers, we must also take into consideration the 
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elementary findings of stability theory (cf., e.g., 10,197). 


It is even possible to use this approximation theory to in- 
vestigate the case of continuous suction as a boundary case with 
a great niimber of single slits arranged in succession. Still, 
a calculation of this sort has not yet been attempted. 

As long ago as 1943 author had presented an imperfectly 
developed form of this theory for suction through single slits in 
a ZWB Report of the AVA Gdttingen (Nr, 1775)Z~2V7, not generally 
available, and correlated it with measurements. In order to cal- 
culate the turbulent boundary layer before and behind the suction 
site, the then-common procedure of E, GRUSCHV/ITZ ^~227 was used. 
It has been shown that, for the correlative measurements in the 
region of Reynold’s numbers that was employed, this procedure 
is only slightly less exact than the newer procedures, such as 
those of E, TRUCKENBRODT and A, WALZ re- 

port we shall again present the essential results of the previ- 
ous calculations and correlations with measurement, but shall 
give only the improved version of the theory, V/ork is now being 
done to confirm the former results with this improved theory, as 
well as to calculate the case of suction of laminar boundary la- 
yers. 


III. 2, A More Precise Statement of the Problem . 

We shall consider a two-dimensional, stationary, incompres- 
sible flow along a solid surface, e.g,, along the contour of an 
airfoil profile. Let us assume that an infinitely long suction 
slit runs transverse to the principal direction of the flow. 

The breadth of this slit is taken to be negligibly small in com- 
parison to the depth 1 of the airfoil, or 



If, now, a lower pressure is used to suction flow mater- 
ial into the interior of the airfoil, this suction flow has 
two diverse effects: 

1 . a direct effect on the boundary layer flow, since 
its material is disappearing into the suction sink; 

2, an indirect effect on the boundary layer, since a 

sink flow is superimposed on the potential flow outside the /12 

boundary layer. The velocity distribution at the edge of the 
boundary layer, altered by the action of the sink, causes the 
boundary layer before and behind the suction slit to develop 
differently, than if the sink were not there. 

A theory for suction must thus comprise, in addition to the 
main problem of boundary layer theory, also a pure problem of 
potential theory, due to the sink effect. In the boundary layer 
calculations, the potential theory distribution of velocity at 
the edge of the boundary layer must be known in advance. For 
this reason we choose to deal at first with this supplementary j 
potential theory problem of the sink effect. It will be seen 
that this effect, even for very small suctions, is considerable. 


III. 3. The Influence of the Suction Sink on the External Velo - 
city Distribution . 

We shall now calculate the additional velocity that a two- 
dimensional, pointlike sink (with linear extension along the air- 
foil span) produces on the surface of an airfoil profile. For 
simplicity we now assume that the profile depth is infinitely 
great and, consequently, that the airfoil surface, upon which 
the sink is situated, is infinitely extended. If Q is the suc- 
tion quantity that disappears in a unit of time across the wing- 
spread b of the airfoil and into the sink, then the potential 

theory additional velocity A u. at a distance a from the sink 

0 
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midpoint is, in the familiar manner. 


A Ug = 


a 

JTab 


( 2 ) 


where Sabu is a cylinder surface with the sink line as axis. 


It is helpful to make the additional velocity Aug nondimen- 
sional with the undisturbed velocity of the approach flow uoo, 
and also to introduce in eq. (2) an airfoil surface (now taken 
to be finite): 


F = b( 


(3) 


where 1 is the depth of the profile, V7e now obtain: 

AUs _ Q — 9 — — = C- -7=J — 

Uqo Uoolb TTa Uoo F Q 5T a 

V/e designate 

a 

''q “ Uco F 


(4) 

m 

(5) 


as the "suction figure". Thus the suction quantity Q is re- 
lated to an imaginary quantity, u«.P, supposed to occur during 
flow with a velocity Uoo along the airfoil surface P in unit time. 


Relation (i|.) can also be used for a finite breadth s of the 
suction slit, if the midpoint of the sink is placed in the mid- 
dle of the suction slit. The calculation in this case proceeds 
only to the edge of the slit, so that the singularity A -*■ oo 
for a^O at the sink midpoint causes no trouble in practice, 

P, EHLERS Z”*2^ (AVA Gottingen) has proved, in an unpub- 
lished work of 1 944> that the simple relation (I4.) for the ad- 
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ditional velocity, resulting from the effect of the siruc, must 
be corrected if a finite wing depth 1 is assumed. In this case 
we must assume a sink on a circular periphery and, in order to 
determine the additional velocity on a given profile, and in- 
clude the sink in the conformal representation of the circle 
( 2 -plane) on the profile contour (^-plane). According to Eh- 
lers, the additional velocity is more exactly expressed by 


_ Cq 
ITR/l 


1^1 (ctg 

'd§l ^ y 




- ctg 


Ta 

2 


) 


( 6 ) 


than by (ij.). Here, alluding to and ^""267, R is the radius 

of the image circle I ^^1 is the absolute value of the differen- 
tial quotient of the image function ^=f(z) that represents the 
circle ( 2 -plane) on the profile plane (5-plane), and ^ is the 
angle on the image circle, with which the individual profile 
points are co-ordinated: = 0 corresponds to the rear edge of 

the profile, <p= (j)^ to the suction site on the profile. 


The total velocity of potential theory, which is de- 
cisive for the development of boundary layers, is, from a linear 
superposition of the main and the additional flow: 


^oo 


A u 


u 


oo 


( 7 ) 


The theoretical investigations of Ehlers have also pro- /1k 
duced the interesting result that the sink flow also Influences 
the circulation about the profile, i.e,, the lift A. This re- 
sult is physically explained by the fact that the sink flow, as 
well as the main flow about the profile, is subject to the 
Kutt a- Zhukovskiy flow-off condition at the rear edge. Prom the 
side of the suction site that is closer to the rear edge, less 
material flows to the sink, since the sharp rear edge hinders a 
flow from the underside of the profile around to the slit. The 
sink flow thus contributes to the circulation about the airfoil 
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profile 


The lift coefficient c is defined, as is known, by: 

ci 

A 

2 GO 

According to Ehlers, the lift coefficient is enhanced by 

Ac„ =2Cg ctg ^ (8) 

due to the effect of the suction sink on the circulation. This 
effect of a sink flow, especially the accelerated increase in 
c as the suction sink approaches the rear edge (^a "*■ 0 ) , has 
been qualitatively confirmed by the measurements of B. REGEN- 
SCHEIT Z~267. For (J)^ = 0, (8) produces the value' ACq = , 

which is no longer physically logical. Evidently formulas (8) 
and also (6) can only be used if and are substantially 

different from zero. In practice, this difficulty in the use 
of the Ehlers formula does not occur, since there are always 
finite widths s of the suction slit. Thus, in the case of for- 
mula (8), the distance of the sink midpoint from the rear edge is 
always s/2 and the corresponding arc on the image circle 
(due to the large value of |dz/d§J in the neighborhood of the 
rear edge) is always relatively large. The additional velocity 
of eq, (6) is only of interest for boundary layer calcu- 
lations up to the edge of the suction slit, or to the distances 
+s/2 from the sink midpoint. Here also the infinity site is 
avoided, ni 

Pig. 10 shows a schematic representation of streamlines 
in the region of a suction slit. An essential feature of this 
pattern is the backwater line (dot-and-dash) which, for a pro-, 
perly chosen suction quantity, ends at the rear edge of the suc- 
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tion slit. The flow material between this line and the profile 
surface disappears into the suction slit. Fig. 11 shows (quali- 
tatively) the corresponding velocity curve at the edge of the 
boundary layer, including the additional velocity due to the ef- 
fect of the sink. It can be seen that the sink creates an ac- 
celeration of the total velocity Ug upstream and a retardation 
doimstream. For boundary layer calculations, at any rate, only 
the velocity distribution above the solid profile surface in 
front of point I and behind point II in fig. 11 is of interest. 
This is the region that can be used with sufficient exactness in 
formula (6). In the region I,II = s directly above the suction 
slit, there are positive additional velocities in the forward 
half and negative ones in the rear half, so that the integral 
of the additional velocities between I and II is equivalent to 
zero. Thus, for the development of the boundary layer directly 
above the slit, i.e. between I and II, the suction sink is, to 
a first approximation, without meaning. Since there is no wall 
friction above the suction slit, we can ignore the finite (and 
presumably very small, according to eq. (1)) breadth s of the 
suction slit for the boundary layer calculations, and take as a 
basis the velocity curve U^(x) shown schematically in fig, 12, 
with a discontinuity at the suction site. 

This discontinuity presents no difficulties for the calcu- 
lations, It is the case, for boundary layer computations, that 
the influences of the acceleration before the slit and of the 
retardation behind are somewhat compensated in their effect on 
the development of the boundary layer profile between I and II, 
when the velocity discontinuity is assumed. 


Ill.i]-. Alterations of the Boundary Layer Parameters by Suction . 
111 . 4 , 1 , The Fundamental Physical Concept , 
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It is now important to ascertain the direct alteration of 
the velocity profile due to a given suction quantity. We can 
assxune that only a portion of the boundary layer material, as 
a rule less than half, is removed. For this case, the stream- /1 6 
line curve in fig. 10 can also be used, A more vigorous suc- 
tion, according to experimental experience, provides no advan- 
tageous effects, and when the auction is to reduce the resis- 
tance it is only a disadvantage, as we have already discussed. 

Fig. 10 already enables us to evaluate how the velocity 
profile is altered by suction: 

A certain velocity profile arrives at the forward edge of 
the suction slit (point I), The backvxater line divides this 
profile into an inner (near the wall) and an outer region. Be- 
hind the suction slit, the inner region of the velocity profile 
has disappeared. The outer region now forms the initial profile 
for the development of the boundary layer behind the suction 
slit. The velocity profile, "clipped" by suction, glides almost 
without friction along the backwater line above the slit and a- 
cross to the rear edge of the latter. In this way, the veloc- 
ity at the "division surface" is reduced to zero when the back- 
water point reaches the rear edge of the slit II. Thus the more 
full velocity profile II behind the slit is derived from the ve- 
locity profile I in front of the slit. At the same time, the to- 
tal thickness of the boundary layer has been reduced from 6j to 
6jj. The exact shape of the velocity profile II is in any case 
unknown. It must therefore be clarified whether and how the ve- 
locity profile II for a given suction quantity and for a given 
velocity profile I caiv be characterized so unequivocally that 
the boundary layer calculation can be continued behind the suc- 
tion slit. We shall limit ourselves here to a boundary layer 
calculation with an approximation procedure, as has been des- 
cribed, for example, in DVL Reports 81; and 136 . In this pro- 
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cedure, the velocity profile of the boundary layer is estab- 
lished by two characteristic magnitudes: a thickness parameter 

and a shape parameter . Our problem, is solved if it is possible 
to calculate the numerical values of these two parameters for 
the velocity profile II from the numerical values of the para- 
meters of the velocity profile I, for given suction quantities. 


III. 4 . 2. Determination of the Boundary Layer Parameters Al - 
tered by Suction . 

III. 1|. 2.1 . A Brief Summary of the Boundary Layer Calculation 
Procedure for Constant Material Values (Incompressible Flow) . 


The above-mentioned approximation theory for boundary la- 
yer calculation basically works with integral conditions /1 7 

for impulse and energy within the boundary layer. These two in- 
tegral conditions are presented in the form of two ordinary dif- 
ferential equations for two characteristic magnitudes of the 
boundary layer. As already indicated, it is convenient to choose 
a thickness parameter and a shape parameter as the character- 
istic magnitudes, and the velocity profile u(x,y) is then gene- 
rally described by an expression in the form: 


u (x, y) 
Ug (x) 


“ ^ U(x) , 



(9) 


Here 6 (x) is the thickness parameter (the total thickness of the 
boundary layer) and H(x) is the shape parameter. The function f, 
which' is different for a laminar and a turbulent boundary layer, 
is assumed to be known in reference to special exact solutions 
or measurements. The definition of the shape parameter can be 
chosen at will. 
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In this approximation theory with integral conditions, the 
following three integral expressions play a special role in the 
case of incompressible flow (to which we limit ourselves here): 



Displacement Thickness Thickness of Impulse Loss 



Thickness of Energy Loss (10,11,12) 

Prom physical and formal reasons it is expedient to intro- 
duce as the thickness parameter of the boundary layer not the 
"total thickness 6 ", which is difficult to define, but the mag- 
nitudes : 


7 = ^ fRf 1" . I, n 5^6*2 

With Rs--]u 


and to introduce as shape parameter: 


n=1 for a laminar 
n=0.268 for a (''3,14) 
turbulent boundary layer 


H = 


6a 


(15) 


These differential equations (with ' indicating the deri- /18 
vative with regard to x) now appear as: 


2'+ Z F/H) - (H) = 0 


(16) 


H' + H-^F(H)“ ^ =0 


(17) 
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with 

F = (2+n) + (1+n)|^(H) j F, = (1 + n) cx( H) 

1 02 ' * 


(18,19) 






( 20,21 ) 


a(H) = partial function of the local frictional coef- 
ficient, dependent on the shape parameter H 
B(H) = portion of the dissipation function dependent 

on the shape parameter H (for a turbulent boun- 
dary layer, e(H)aB =0.0056) 


laminar boundary layer 

n = N“1 

('22 ) 


n = 0.268, ■ , 


turbulent boundary layer 

1 n-N = 0,1 

N = 0.168 

(23) 


u^ 5 (x) is the velocity distribution, to be assumed as given, 
at the edge of the boundary layer; in the case treated here, it 
must also include the additional velocity of the suction sink. 

In DVL Report 136 analytical expressions for the uni- 

versal functions to for laiuinar and turbulent boundary la- 
yers (for incompressible and compressible flow, with and without 
heat transfer) as well as solution methods for the system of e- 
quations (l6) and (17) have been given, so that no further ex- 
position is needed here, 

V/e must now answer the question as to how the parameters 
Z and H during incompressible flow are altered under the influ- 
ence of an unsteady suction by a slit. It should also be men- 
tioned that there are no fundamental difficulties in generali- 
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zing this theory for compressible flow. 


Derivation of the Connection between the Suction 
Number C„ and the Boundary Layer Par ameter s. 

' y ^ ^ 

We have already established in section III.Ii-.l. that the 
shape of the velocity profile II behind the suction slit is not 
at first completely known. But since, within the framework of 
the approximation theory sketched in section , , inte- 

grals of the form (10), (11), and (12) on the velocity profile 
are essentially concerned, it is possible to represent these 
integral expressions at the site II, i.e. behind the suction, 
by the difference of two integral expressions. StarWng from 
the physical representation of the flow curve in the region of 
the suction slit, as' shown in fig. 11a, the following are valid: 



(24) 


(25) 

( 63)11 "" 

( 26 ) 


with 


y<i 



( 27 ) 



(^)q =/ V('-^)dy 

( 28 ) 

°ya 



(29) 


0 


Here y^ is the distance from the wall to which the veloci 
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ty profile I, arriving at the suction slit (fig. 11a)# must be 
"clipped” in order that the given suction quantity per unit of 
time 


Q= u„bl 


( 30 ) 


disappears into the suction slit. With (30), 

yq 


Q= Cq Uoo bl = 




(31) 


must be valid. Thus the "thickness of suction" 6 q is defined as: 




ya 


/20 

(32) 


Prom (31) and (32) it also follows that 


c>~ g _ - _h2. , 

As we shall later show in detail, there is no difficulty in 
calculating the integral expressions (2?) to (29), and thereby 
the values (<S-|)jj» ^ given 

group of laminar or turbulent velocity profiles according to eq. 
(9) as functions of the suction number c„, the shape parameter 
H, and the local velocity ratio u^/u^o . 

In this way, from eqs, (13) and (lij.), the thickness para- 
meter Zjj can also be immediately determined as the initial va- 
lue of the computation behind the suction slit. 


It appears at first somewhat more difficult to ascertain 
the shape parameter of the velocity profile, altered by 
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suction, from the magnitudes i.e*> to 

again classify this somewhat ‘’truncated" velocity profile in the 
form category of expression (9). 


For a single-parameter family of curves, which we have 
taken as a basis in the approximation theory with the general 
expression (9)> the relationships ^3^ ^2 ~ ^ 

(15)) are specifically correlated. Thus ^ specific 

function of the chosen shape parameter H. For the velocity pro- 
files altered by suction we can no longer assume that they may 
be accurately classified in the original family, V7e must much 
rather expect that the relationships 


Tsir'"“ 

^ II 


(33) 


of the profiles altered by suction can no longer be specifically 
coordinated to each other. If we nonetheless assume the expres- 
sion (9) of the original velocity profiles to be valid, then we /21 
obtain, from (<S-i )u/(<S2^II rather dif- 
ferent HTT-values: H-pt a^id Htt • It is now our concern to 

-LX XXi 2 XX3 2 

arrange, at least approximately, the velocity profile altered by 
suction in the original^ family of single-parameter velocity pro- 
files, i.e., to find in this family a profile that is most simi- 
lar to that altered by suction. A physically reasonable pre- 
scription for this arrangement of the "clipped" velocity pro- 
files in the original family of curves is now, evidently, to form 
an arithmetic mean from the two generally different values 
and H, 


Hll 


12 


^TT * 
J.132 


"n = 


H„ + 


H 


II 


32 


(3i^.) 


and to continue the calculation behind the suction site II with 
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This principle for the determination of is the 


this value, 
more dependable as the values of Hjj 
ferent from each other. 


and H 


12 


II 


32 


are less dif- 


As we shall later establish, this prerequisite for the de- 
pendability of the approximation theory is fulfilled in a use- 
ful manner. 


IV. The Application of Suction Theory to Turbulent Boundary 
Layers . 

IV. 1 . Usage Formulas on the Basis of Power Profiles . 

We shall now occupy ourselves at first with the suction of . 
turbulent boundary layers. For the velocity profiles in this 
case (at any rate for the calculation of the integral expres- 
sions ( 24 ) to (29))# we can take as a basis the simple exponen- 
tial expression 


( 35 ) 

with 

0<k<1 or more exact 0.1 <k<0.7 (36) 

On this basis, the integral expressions (10) to (12) can be im- /22 
mediately and continuously evaluated. These are the expressions 
that are valid for the velocity profile in front of the suction 
site, and which we must thus distinguish with the index I. 

From (10), (11), (12) with ( 35 ) it follows that 


22 



I . 

^ 6 l+kj ^ (1+ kj)(l + 2ki) ' 


) = —I 

^ 5 (1 + kj )(i +3kj.) 


(37,38,39) 


From which is obtained 


^2 


5-1- 


iko) 


and, with (1^), 
6 


(— ) = 

^ 62 ^ 


(6z) 

Si 

(5a) 


, 1*2ki 2 -Hj 

^ l+3kj kj- 3 h^-a 

ikUk2) 

ons ( 24 ) to 

( 29 ) assume the form. 

with (35), 


l+k_ 



(yfl/<5i ) 


0)1 
1— 1 

1 

1 + kj. 

0^3) 

(yaMi )'*''! 

(yQ/6j- ) 

m) 

1 + kj 

1 + 2kj 



0l5) 


1 + k. 


l+3k. 


From this it follows, for the magnitudes behind the suction site 
(distinguished by the index II), that 


(S,)n 


— 


1+k, 


1+kj 

Va ^ ^ya/6i) 

6t 1+k 


01-6) 
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( 47 ) 



_ 


(Yq/^i) 

5i 

(1+ki)(1+2kj) 

1 +k 
I 

1 + 2k 

I 

(63 )„ 

_ 2 k, 

(yoMi) 

1 + 3kT 

(ya/5i) 

5i 

“ (1 + ki)(U3l<j) 

1 + kj 

l+3kj 


(48) 


For the thickness of suction 
that 


with ( 32 ), 


it is now valid 




1+k, 


(yn/<si) 


1 + k. 


(49) 


From which the magnitudes, needed in (46) to (48), with (32a) 
and ( 38 ) yield 






(<5; / 6 )n 

( J 


1 + k- 


r -I 

L (u^/u^)^ H 2 kj (62/Dj J 

All the magnitudes appearing in (50) are known: kj(Hj) fol- 
lows from (4l)j I according to eqs, ( 13 ) and (l4)j 

and u6/uoo and Cq are given in the original statement of the prob- 
lem. 


We can now, to check the methods, described in section 
III.4«2.2, for the determination of the shape parameter be- 
hind the suction slit, easily form the relations (^-)/^2^II 
( 62 / 62 ^ II* result is shown in fig. 13* with kj(Hj) as 

parameter is diagrammed above the "suction distance" 7 q/^j 
( which is connected to the suction niunber c^, according to ( 50 ))» 

On the one hand, with (47) and (48), follows directly from /24 
the definition (33): 
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( ^3 )ll 


= H 


II 


32 


On the other hand, one can calculate Hjj from t>y way 

of relations (J 4 .O) and (i |.1 ), namely as: 



1+ 3kjj 


(51) 


with 


( 51 a) 

As can be seen from fig. 13» the values H-- and 11^^- are 

only slightly different, so that we should regard the value 
obtained by arithmetic averaging according to eq. ( 2 i|.), as a 
usable initial value for the shape parameter behind the suction 
site. 


This analysis is immediately valid for turbulent boundary 
layers. For laminar boundary layers a corresponding verifi- 
cation is still wanting. 

The value which in addition to is used to continue 

the boundary layer calculation behind the suction site, follovjs 
unequivocally from eqs. (]+ 7 )> ( 13 )» and ( 14 ). 


IV. 2. A General Appraisal of the Suction Effect . 

With the simple derived relationships, we can now make a 
quick survey of the effect of a definite suction quantity on the 
boundary layer. For this purpose, we begin with fig. 1 3» from 
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v/hich it is immediately apparent to what relative distance from 
the wall, y^/6-, the velocity profile with shape parameter 
Hj(kj), arriving at the suction slit, must be "clipped" in or- 
der for a velocity profile with a prescribed value Hjj to arise 
behind the suction slit. The relative suction distance 7^/^^ 
and the suction nxunber c^ are connected to each other by eq, ( 50 ). 
The other magnitudes appearing in eq, (50)--ug/uoo, k(Hj), and 
( 62 / 1 ) 2 — are known with the calculation of the potential flow 
and the boundary layer up to the suction site I, The relations 
(i|1 ) and ( 42 ) provide the connection between k and H, /25 

We shall discuss the application of eq, (^0), in combina- 
tion with fig, 13> by an example. For this example vre shall 
assiime that the turbulent boundary layer in front of the suc- 
tion slit has passed through a region of pressure rise of such 
intensity that, without suction at this location, flow detach- 
ment would immediately ensue. According to detachment of 

the turbulent boundary layer can be expected when the shape pa- 

*1 

rameter H is somewhat less than the value 1.57 : 

Hflet. "1-57 (52) 

To this value of there corresponds, according to eq, (42), 

the value 


'=det. " <53) 

We should now deteimiine the suction quantity that reduces 
the shape parameter Hj = 1*57 (k^ = 0 , 606 ) to the value Hjj = 


-I 

In many instances flow detachment is first observed when H < 

1 ,53» Such cases apparently prevail when a pressure drop, fol- 
lowing a pressure rise, suctions away the dead water region, 
formed by the detachment, in the flow direction. 
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1.80 (kjj = 01,43). This value of characterizes a veloci- 
ty profile such as that appearing in a flow without pressure 
gradients and with large Reynold’s n\ambers (cf. J[7> DVL Re- 
port 136, p, 51 )• To achieve such a suction effect, one must 
"clip*' the velocity profile to Yq/ = 0,53» according to fig. 
13. The appropriate suction number c^ is, from eq, (50): 



I + k 




6 * 1.606 
-^) = 0,53 


2.212 

0;606 I 1 


= 1.32 ( 




I 


) 


I 


(54) 


The suction quantity Q then follows from eq, (5). 


Thus Cq is proportional to the thickness of impulse loss 
62/1 (related to the wing depth 1), For large Reynold’s num- 
bers, as can be assiuned in the region of validity of boundary 
layer theory, 62/I is a very small niimber (almost one power of /26 
ten smaller than the ratio 6/1). Thus, c„ is also very small. 

For Reynold’s niombers R-, = pu«l/y of magnitude 10 , c„ in the 

chosen example has the magnitude of 10 , according to calcu- 

2 

lations and measurements. 


We shall now evaluate the sink effect for c^ = 0,01, using 
the simple formula (4). For c^ = 0,01, the following additional 
velocities AuyUoo obtain, dependent on the distance a/1 from 
the sink midpoint on the profile surface: 


2 

It should be mentioned here that, when using suction to main- 
tain the laminar! ty of the boundary layer, even smaller values 
of c„ come into consideration. In the experimental findings of 

Pfenninger [_ 1 67. 9 cited at the close of section II, c„ has the 
magnitude of 1 O”^, ^ 
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a/I 

0 

0.01 

0.0 2 

0.05 

0.10 

^ > 
oo 

oo 

0.318 

0.1 59 

0.0 6 A 

0,032 


It is seen that the additional velocities at = 0.01 are al- 
ready so large that they must be taken into account in the cal- 
culation of the boundary layer. 


V. ^ Example from Aerodynamics ; Ascertainment of the Host 
Favorable Disposition of the Suction Slit for Increases of 
°a max — NACA - profile 2301^ . 

If it is intended to suppress flow detachment by suction, 
then the question arises as to the disposition of the suction 
slit so that the suction quantity for a given velocity curve 
u6(x) becomes a minimum. 

To begin, we can make the following rough predictions a- 
bout the slit disposition that can bring ,about a minimum suction 
quantity; 

The suction slit should be disposed downstream at the lat- 
est where flow detachment would ensue without suction (with uti- 
lization of the sink effect, which somewhat reduces the pressure 
rise, also slightly behind this position). But the slit should 
not lie to much in front of the detachment site, since a boun- 
dary layer, made thinner by suction, grows more quickly than a 
thick one, due to the increased friction. /27 

The most expedient disposition of the suction slit can now 
be clarified by the developed theory in combination with a cal- 
culation for pressure distribution and boundary layer. 
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As an example for such a calculation, in 1942 in the Aero- 
dynamical Research Institute, Gottingen, an airfoil, profile of 
the type NACA 23015 was chosen (ZV/B Report 1775) Z~2l7. To 
check the theory, in a water channel of the AVA Gottingen mea- 
surements of the pressure distribution were carried out for this 
profile, along with flow observations by the spread of a dye. 
From the measurement of the pressure distribution p(x)/ ~ , 

the velocity distribution u6(x)/u<„ , necessary as a basis for 
the boundary layer calculation, was obtained from Bernouilli ‘ s 
Principle. The flow observations enabled recognition of both 
the turning point and the detachment point of the turbulent 
boundary layer, with good accuracy. In this way, the prerequi- 
sites needed to compare measurements with boundary layer calcu- 
lations (with suction) were produced. 


At that time the integral condition of the impulse, eq. 
(l6), for laminar and turbulent boundary layers, was already 

A 

available for the boundary layer calculation. The "wall bind- 
ing" served as the second equation for a laminar boundary layer 
(see, e,g,, Z”*21.7, eq. (11a)). For a turbulent boundary layer, 
as already mentioned at the beginning, the Gruschwitz differen- 
tial equation for the shape parameter was used Z~2l7 (of. eq, 

( 38 )): 



i!_)’ =) _ (Aa. 


2k . 

) = 1 - [• 


, 2h 


( 1+k) (1 + 2k ) 


( 55 ) 


1 

For a turbulent boundary layer, the empirical law for the shear- 
ing stress of the wall was not yet so accurately known as is to- 
day possible, after the works, e.g,, of Ludwieg-Tillmann Z"~297. 

In calculating the detachment point of turbulent boundary la- 
yers, this uncertainty was not really significant, since in this 
case the frictional member is subordinate to the pressure member 
in the pressure equation. 
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Eq, (42) provides the connection with the magnitude H, used by 
us. K. WIEGHARDT ^~"2j7 has proved that this equation, with con- 
stants that have been empirically determined from a limited num- 
ber of measurements, can be interpreted as an integral condi- 
tion for the energy. The equation has shown itself to be very /2d 
useful in the region of Reynold's niunbers for which the con- 
stants had been determined. As a rule, its results are only 
slightly worse than the new energy equation (17)» which has a 
better physical foundation (cf., e.g., J. ACKERET ^ 287) • 

Thus the comparison between calculation and measurement, 
carried out in the cited internal report, is essentially still 
valid. We therefore reproduce several pictures from this re- 
port, which reveal the essential findings (figs. 1i|. to 19). It 
can be established that the boundary layer calculation, with con- 
sideration of the suction influence, renders rather well the 
practically attainable c. values as a function of the sue- 
tion nixmber c^, and the angle of attack a. The theory also pre- 
diets well the most favorable disposition of the suction slit 
(in the present profile example, x/1 = 0.2), as seen in figs. 21 
and 22. 

A check with the somewhat improved theory, presented in this 
report, is in preparation. The improvements, as already men- 
tioned, concern the integral condition for the energy, eq. (17)> 
and furthermore the determination of the shape parameter Hjj 
behind the slit. In the earlier calculations, only the relation 
^^1 ^II'^^ ^2^11 used to determine the shape parameter of that 
time, q(H). According to the improvement of the theory pre- 
sented here, the shape parameter (H)jj of the velocity profile 
can be determined more certainly as a mean from the relations 

(61)11/(62^11 (34)* As seen in 

fig. 13, this mean H does not differ essentially from its two 
components. This fact also indicates that nothing decisive will 
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be altered when the former result of the comparison between 
measurement and calculation is subjected to verification. 

The result of such a verification, as well as the application 
of the here-depicted theory in the case of a laminar boundary 
layer (the problem of maintaining the laminarity of the boun- 
dary layer by numerous, successive, single slits; see, e.g., 

the measurements of Pfenninger shall be reported in 

2 

due time. 


VI, Summary . /29 

The majority of the presently known theories for the suction 
of the flow’ s boundary layer concern the physically optimal case 
of the continuous distribution of the suction quantity over the 
surface, which is at the same time mathematically more easy and 
more exact. In practice, at any rate, to approximate the phy- 
sically optimal case of the suction effect, reasons of design 
dictate a rather large number of successive single slits, but 
also and frequently only a single slit, if, e.g., flow detach- 
ment is to be avoided. An exact treatment of this case is dif- 
ficult, On the contrary, an approximate solution, especially 
in connection with the calculation procedures for laminar and 
turbulent boundary layers, described in DVL Reports 8i|.and I 36 , 
can be relatively easily produced. The author had already de- 
veloped an imperfect form of this suction theory in 19ll3» Tbe 
most important results of this report, at the time enjoying only 

2 

During the printing of this report, VJ. V/uest, AVA Gottingen, 
has programmed the approximation theory, developed here, for a 
turbulent boundary layer. Systematic calculations with many 
successive single slits have produced a satisfactory accord 
with experiment. 
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limited publication, have been again presented in the current 
report. Comparisons with measurements have demonstrated the 
fundamental validity of this simple theory, above all for turbu- 
lent boundary layers. Today this theory can also be checked 
in the case of suction to maintain the laminarity of the boun- 
dary layer, in combination with the stability theory and the 
measurements of W. PPEMINGER, Works in this direction are 
being prepared. 

For suction by single slits, the external velocity field 
of potential theory, i,e, the starting basis for the boundary 
layer calculation, is influenced by the suction sink. Simple 
evaluation formulas are adduced for this influence, 

Alfred Walz, University Lecturer and 
Doctor of Engineering, Freiburg, FRG. 
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Figures 




(Theohe SchLichting) 


Figures 3 and 4 * Key; a) suction; b) suction slit; c) con- 
tinuous suction (Schlichting Theory) 
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Figure 5. Key: a) blowing 



Figure 6, Ordinary flow 
in an abruptly broadened 
channel. 



Figure 7, Flow in an ab- 
ruptly broadened channel 
with suction at the walls. 
The white markers indicate 
the sites of the (invisi- 
ble) suction slits. 











Figure 9 . Key: a) 

entj b) even platej 


total resistance coeffici- 
c) transition lara-turb. 
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Figures 10, 11, and 12 
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Figure ll|. 
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Figure 16. 



21° i Cq » 0,0050. 


Figure 17 
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Figure 20, Key: a) slit; b) without suction; 

c) (calculation for even flow);^ d) (measure- 

ment); e) measurement (PB Nr, 1611); f) calculation; 
g) geometrical angle of attack; h) effective angle 
of attack. 
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